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摘 要 : 干旱 灾害 是 京 津 沉 地 区 主要 气象 灾害 之 一 ,严重 影响 社会 经 济 和 粮食 安全 。 评 估 京 津 浊 地 


区 干旱 灾害 危险 性 ,识别 高 危险 区 ,站 在 为 灾害 防范 和 应 急 调控 提供 参考 。 利 用 1958—2017 年 37 


个 气象 站 降水 日 数据 ,计算 标准 化 降水 指数 (SP16) ,根据 游程 理论 界定 干旱 历时 和 干旱 强度 。 运 
用 K-S 检验 法 AIC 法 分 别 确定 指标 的 最 优 边 缘分 布 函 数 和 Copula 函数 ,计算 5a10a、30a 及 50&a 
一 过 水 平 下 干旱 强度 和 干旱 历时 的 重 现 期 。 结 果 表 明 : 京 津 费 地 区 20 世纪 60 年 代 中 期 至 今 SP16 
旦 下 降 趋势 ;指数 分 布 和 广义 极 值 分 布 分 别 是 干旱 历时 和 干旱 强度 的 最 优 边 缘分 布 ,Frank Copula 
在 变量 联合 概率 分 布 计算 时 拟 合 最 好 ;就 干旱 历时 和 干旱 强度 重 现 期 而 言 ,30 a、50 a 一 遇 水 平 下 
联合 重 现 期 为 ,北京 市 .沧州 市 .衡水 市 及 承德 市 北部 等 地 区 的 重 现 期 长 ,危险 性 低 ,其 余地 区 危险 
性 较 高 ;站 点 干旱 灾害 的 联合 重 现 期 理论 值 与 实际 情况 较 符合 ,能 够 反映 实际 灾害 情况 。 
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全 球 变 暖 背景 下 ,极端 气候 事件 ,尤其 是 近 几 十 
年 来 的 干旱 ,高温 等 …… 极端 气候 灾害 频 发 ” 。 干 旱灾 
害 作为 严重 影响 全 球 的 自然 灾害 之 一 ,危害 范围 大 ， 
持续 时 间 长 。 到 21 世纪 末 , 全 球 干旱 半 干 旱 区 将 占 
陆 表 面积 的 30% 以 上 ,其 中 75% 将 发 生 在 发 展 中 国 
家 。 中 国 每 年 因 干 旱灾 害 造 成 的 损失 巨大 ,民政 
部 、 国 家 减灾 委 统 计 ,2017 年 造成 农作物 受灾 98 748 
km? ,绝收 7 524 km? ,直接 经 济 损 失 3.75 x 10" 元 。 

中 国 干旱 灾害 多 发 地 主要 分 布 在 华北 . 黄 淮 \ 西 
北 东 部 .东北 西部 .华南 西部 及 西南 大 部 等 。 而 
位 于 华北 平原 的 京 津 费 地 区 ,干旱 灾害 形势 严峻 , 干 
早 覆盖 率 除 冬季 外 均 呈 不 同 程度 的 增长 趋势 ”1 。 
降水 呈现 轻微 减少 趋势 ,而 温度 升 高 趋势 明显 ,这 共 
同 作用 导致 了 明显 的 干旱 化 “1。 干旱 灾害 不 仅 影 
响 社 会 生产 生活 ,还 会 影响 农作物 生产 ,威胁 河北 省 
力 至 全 国 的 粮食 安全 '" ,这 对 京 津 蔓 地 区 干旱 灾害 
防 灾 减 灾 提出 了 新 的 要 求 和 挑战 。 

现 有 研究 中 ,干旱 灾害 识别 主要 依托 干旱 指标 
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进行 界定 。 由 于 致 灾 机 理 复杂 ,在 界定 时 不 仅 要 考 
虑 降水 ,还 要 考虑 蒸 散 、 土 壤 含 水 .温度 及 径流 等 。 
目前 ,研究 和 应 用 的 气象 干旱 指标 众多 ,包括 降水 距 
PHSR .相对 湿润 度 指 数 .Z 指数 .标准 化 降水 
指数 ”~ 中 土壤 相对 湿度 干旱 指数 、 帕 默 尔 干 旱 指 
数 、 连 续 无 有 效 降水 日 数 有 效 降水 指数 .综合 气象 
干旱 指数 及 标准 化 降水 蒸 散 指数 哺 ] 等 。 其 中 ， 
标准 化 降水 指数 (SP1) 计算 简单 ,干旱 监测 与 预测 
能 力 较 强 ,对 干旱 变化 反应 敏感 ,能 够 较 好 地 反映 干 
旱 历时 和 干旱 强度 ,被 广泛 应 用 于 干旱 监测 与 评估 
中 中- 由。 灾害 事件 的 发 生 是 多 变量 相互 联系 共同 
作用 的 结果 ,因而 在 灾害 风险 研究 中 ,多 维 致 灾 因 子 
评估 逐渐 发 展 和 应 用 。 目 前 多 维 联合 分 布 的 方法 研 
究 中 ,Copula 函数 由 于 其 对 变量 的 边缘 分 布 无 要 求 ， 
变量 可 存在 相关 关系 ,灵活 性 和 适用 性 较 强 ,广泛 应 
用 于 水 文 干旱 灾害 研究 中 ~" 。 在 干旱 灾害 研究 
中 ,以 干旱 历时 和 和 干旱 强度 的 二 维 联合 研究 最 多 , 近 
年 也 逐渐 发 展 到 三 变量 REKET 。 
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其 中 ,中 国 研究 区 主要 涉及 西南 地 区 '" 、 黄 土 高 
原 " 珠江 流域 "中 汉江 流域 "及 广东 省 中 等 。 

干旱 灾害 特别 是 重大 干旱 灾害 监测 和 防御 是 科 
学 研究 的 优先 主题 。 科 学 地 预测 干旱 灾害 趋势 , 提 
出 行 之 有 效 的 对 策 , 以 适应 国家 防 灾 减 灾 、 应 对 气候 
变化 和 保障 粮食 安全 的 迫切 需求 。 京 津 晕 地 区 作为 
干旱 重 灾区 之 一 ,其 干旱 灾害 时 空 变化 及 成 因 已 为 
诸多 学 者 研究 5 ,但 尚未 有 针对 该 地 区 利用 Copula 
函数 开展 干旱 历时 和 干旱 强度 二 维 概率 及 重 现 期 的 
相关 研究 。 因 而 本 文 基于 标准 化 降水 指数 ,借助 
Copula 函数 探究 京 津 蓉 地 区 干旱 灾害 历时 和 强度 的 
联合 分 布 ,分 析 其 时 空 特征 ,识别 干旱 灾害 高 危险 
区 ,以 期 为 该 地 区 干旱 灾害 防 灾 减 灾 能 力 提 升 .区 域 
经 济 发 展 及 粮食 安全 保障 提供 参考 。 


1 研究 区 概况 及 数据 方法 


1.1 研究 区 概况 和 研究 数据 

京 津 费 地 区 位 于 113°04’ ~ 119°53'E,36°01' ~ 
42°37'N, 具 有 四 季 分 明 、 雨 量 集中 日 照 充 足 及 干 湿 
期 明显 等 特点 。2017 FREAK AAH 1.11 x 
10* 人 ,农作物 总 播种 面积 93 825 km’, GDP 合计 
8.255 9 x10 元 , 占 全 国 的 9.98% ,是 中 国 经 济 的 
第 三 增长 极 , 同 时 也 是 国家 政治 文化 中 心 和 北方 经 
济 中 心 。 

考虑 长 时 序 气象 数据 的 完整 性 ,确定 研究 时 段 
为 1958—2017 年 ,涵盖 了 京 津 焉 地 区 37 个 国家 气 
象 站 (图 1) 。 数 据 类 型 为 降水 逐日 数据 ,来 源 于 河 
北 省 气象 局 和 中 国 气象 数据 网 (http: // data. cma. 
cn/) 。 粮 食 产 量 和 面积 数据 来 源 于 《河北 经 济 年 鉴 
2016 兴 中 国 统计 年 鉴 2016》。 
1.2 研究 方法 

人 研究 方法 主要 包括 SPI 计算 、Mann Kendall 突 
变 分 析 ,游程 理论 及 Copula 函数 联合 分 布 计 算 。 不 
同时 间 尺 度 的 SP7 反映 不 同 尺 度 的 干旱 灾害 情况 ， 
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指数 ( SPI) 事件 识别 


图 1 国家 气象 站 点 分 布 


Fig.1 Distribution of national meteorological stations 


1 个 月 尺度 的 SPN 能 较 好 地 反映 降水 引起 的 旱 涉 
频繁 交替 ,3 个 月 和 6 个 月 尺度 的 SPB .SP16 对 旱 涝 
暴发 与 持续 时 间 的 辨识 更 敏感 ,而 12 个 月 尺度 的 
SP112 对 长 期 旱 涝 影 响 具有 显著 的 标识 功能 ”| 。 为 
了 更 好 地 反映 干旱 对 京 津 蔓 地 区 的 长 期 影响 ,选取 
对 干旱 持续 时 间 辨 识 敏 感 , 且 可 以 反馈 长 期 干旱 的 
SP16。 王 旱灾 害 重 心 的 时 空 分 析 中 ,考虑 干旱 灾害 
是 一 个 长 期 存在 ,缓慢 演变 的 自然 灾害 ,将 研究 时 段 
(1958 一 2017 年 ) ,结合 我 国 社会 经 济 发 展 的 重要 时 
期 ,以 20 a 为 基准 ,划分 为 改革 开放 前 20 a(1958 一 
1977 年 ) 改革 开放 后 20 a(1978 一 1997 年 ) .进入 
新 世纪 以 来 的 近 20 a( 1998—2017 年 ) 。 借 助 游程 
理论 识别 干旱 历时 和 和 干旱 强度 ,分别 确 定 其 最 优 边 
缘分 布 ,代入 Copula 函数 计算 重 现 期 ,基于 ArcGIS 
实现 变化 特征 的 空间 表达 ,评估 京 津 蔓 地 区 干旱 危 
险 性 (图 2)。 


Copula 时 空 干旱 灾害 危 
联合 分 布 分 布 险 性 评估 


图 2 京 津 豆 地 区 干旱 灾害 危险 性 评估 技术 路 线 


Fig.2 Technical framework of drought hazard assessment in the Beijing-Tianjin-Hebei region 
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1.2.1 干旱 灾害 识别 及 指标 界定 ”借助 游程 理 
论 趾 确定 干旱 灾害 事件 及 其 特征 指标 ,SP16 参照 国 
家 标准 《气象 干旱 等 级 GB/T20481 一 2006》 进 行 计 
算 及 等 级 划分 ( 表 1)。 


表 1 不 同 干旱 类 型 的 SPI 取 值 
Tab.1 Drought types corresponding to the SPI value 


等 级 类 型 SPI fi 

1 无 旱 -0.5 < SPI 

2 轻 旱 -1.0<SPI< -0.5 
3 中 旱 -1.5 <SPI< -1.0 
4 重 旱 -2.0<SPI< -1.5 
5 特 旱 SPI< -2.0 


一 次 干旱 灾害 事件 涉及 干旱 历时 ( 忆 ) 和 干旱 强 
度 (5)。 干旱 历时 为 一 次 干旱 事件 的 持续 时 间 , 干 
晶 强度 是 干旱 历时 内 SPI 的 累加 , 即 图 3 中 阴影 音 
分 的 面积 ”2 。 分 别 选 取 SPL BYR R = -1\R, = 
-0.5 R, =0 ,干旱 特征 指标 界定 的 过 程 为 : 

(1) 若 SP1 小 于 届 ,初步 判断 此 月 为 干旱 , 即 
图 3 PIN Be a,b de; 

(2) 当 干 旱 历 时 为 1 个 月 (ae) 且 SP 小 于 
Ry WIA 1 次 干旱 事件 (e) ,反之 不 计 为 干旱 事件 
(a); 

(3) 若 相 邻 2 次 干旱 事件 间隔 为 1 个 月 (b 和 
d) , 旦 间隔 月 份 的 SP/ 小 于 R,, 即 S$. <0, 则 合并 为 1 
次 干旱 事件 ,干旱 历时 D =D, +D, +1, 干 旱 强 度 
(阴影 面积 )$ =S, +S, ,否则 为 2 次 干旱 事件 。 
1.2.2 按 缘 分 布 函 数 构 建 ” 首先 确定 干旱 历时 和 
干旱 强度 的 边缘 分 布 ,选用 水 文 和 气象 领域 常见 的 
9 种 单 变量 分 布 孔 数 , 即 伽 马 、 指 数 、 对 数 正 态 、 威 布 
AR ERS 广义 极 值 . 极 值 . 泊 松 及 广义 帕 累 托 分 布 。 
利用 Kolmogorov -Smirnov 检 验 确 定 最 优 边缘 分 布 。 


d 干旱 强度 ° 


到 3 游程 理论 识别 干旱 特征 变量 


Fig.3 Identification of drought characteristics using run theory 


同时 采取 Kendall 秩 相关 系数 和 Spearman 秩 相 关系 
数 度量 变量 间 的 相关 性 ,参数 估计 采用 极 大 似 然 法 。 
1.2.3 Copula 函数 选择 及 其 参数 估计 Copula 方 
法 最 早 于 1959 年 提出 ,理论 基石 是 著名 的 Sklar 定 
理 ,可 以 表述 为 : 设 X,Y 为 连续 的 随机 变量 ,其 边缘 
分 布 函 数 分 别 为 Fy (x0) A Py Cy) ,联合 分 布 函数 为 
F(x,y) ,大 各 边缘 分 布 函 数 连续 , 则 存在 唯一 的 函 
数 Copula 函数 C (u,v) ,使 得 

F(x,y) =CoCEFCxz) rr)，Vzy (1) 
式 中 :f(x,y) 是 联合 分 布 函数 ;Cs 是 Copula 函数 ;6 
为 待定 参数 。 

Copula 函数 大 体 上 可 分 为 椭圆 型 Archimedean 
型 及 二 次 型 3 种 类 型 ,其 中 Archimedean Copula 以 
其 构造 简单 . 易 求 解 等 特点 应 用 广泛 。 目 前 干旱 分 
析 中 常用 的 Archimedean Copula ,包括 Clayton 族 、 
Frank 族 及 Gumbel 族 ( 表 2)。 


表 2 三 种 Archimedean Copula 函数 表达 式 
Tab.2 Three bivariate Archimedean Copula functions 


Copula 


参数 取 值 范 围 


Clayton max[ (utei) 0 de[-1,0)\{0} 


-0u _ = 
Frank 一 Pi ie dER\{0} 


ety 


Gumbel opf - [Cm i) |" de[1,~) 


采用 均 方 根 误差 (RMSE) 和 AIC 信息 准则 确定 
最 优 Copula 函数 。 
1 


MSE = ~~ >, (Pe, - P) (2) 
AIC =nlog /MSE +21 (3) 
RMSE = \/MSE (4) 


式 中 :Pei 为 干旱 联合 经 验 概率 ;P; 为 Copula 函数 联 
合 概率 ;n 为 样本 数 ;i=1,2…,n;l 为 模型 中 参数 个 
1.2.4 和 干旱 指标 重 现 期 确定 ”干旱 历时 ( 刀 ) AMF 
旱 强 度 ($) 的 单 变 量 重 现 期 为 : 

M 


Py = Te (gy ho (5) 
M, 
Ts=T IF sys) (6) 


式 中 :Ff,(d) 和 了;(s) 为 干旱 历时 和 干旱 强度 的 边 
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FRAY Api KRM, 为 2 次 干旱 事件 的 平均 时 间 。 
E(L) 


E(L) 


二 维 联 合 重 现 期 与 同 现 重 现 期 分 别 为 : 
E(L) 


T(s,d) = 


E(L) E(L) 


T'(s,d) = 


P(SSsUD2d) [1-F(s,d)] [1-C(Fs(s),F,(d)) ] 


(7) 


E(L) 


式 中 :7(s,d) 为 联合 重 现 期 ;让 (s,d) 为 同 现 重 现 
期 ;E(L) 为 连续 事件 的 间隔 期 ;fp(d) 和 Fs(s) 分 别 
为 干旱 历时 和 干旱 强度 的 边缘 分 布 函 数 。 


2 结果 与 分 析 


京 津 责 干旱 年 际 变化 特征 
1958—2017 年 , 京 津 囊 地 区 站 点 SP16 平均 值 
呈 上 下 波动 (图 4) , 共 发 生 干 旱 166 次 (SPI<0.5)。 
对 于 时 间 尺 度 较 长 的 SP1, 多 次 持续 降水 才 会 影响 
其 波动 1]。 借 助 Mann-Kendall 进行 突变 分 析 。 当 
UF >0 ,表明 序列 呈 上 升 趋势 ,UF <0 表明 序列 呈 下 
降 趋 势 , 超 过 临界 值 ( +1. 96 ) 表 明 趋 势 明 显 。 如 果 
UF 和 UB 曲线 出 现 交 点 且 在 临界 值 之 间 ,那么 交点 
对 应 的 时 间 就 是 突变 时 间 。 图 5 中 ,20 世纪 50 年 
RACES HX SP16 呈 上 升 趋势 ,20 世纪 60 年 代 
初 呈 下 降 趋 势 ,1964 一 1966 年 又 转 为 上 升 , 在 20 tH 


2.1 


P(S>snD>d) [1-Fy(s) —F)(d) +F(s,d)] [1 -Fy(s) —F,(d) +C(Fs(s) ,F,(d)) ] 


(8) 


纪 50 年 代 末 ~ 60 年 代 初 出 现 多 次 突变 。1966 年 至 
今 整体 下 降 , 其 中 1979—1980 年 短暂 上 升 。 在 20 
世纪 80 年 代 、1992 一 1994 年 .90 年 代 末 ~21 世纪 
前 10 a 下降 最 显著 ,2015 年 左右 又 出 现 一 次 突变 ， 
变化 趋势 与 王 宏 等 *| 的 结论 基本 一 致 。 

在 1958 一 1977 年 .1978 一 1997 年 及 1998—2017 
年 这 3 个 时 期 中 ,1978 一 1997 年 间 发 生 干 早 频 次 最 
高 为 60 次 ,其 中 轻 度 灾害 最 多 为 45 次 。 而 1958 一 
1977 年 间 发 生 干 旱灾 害 51 次 ,1998 年 后 20 a, 发 生 
干旱 灾害 55 次 ,这 两 个 阶段 中 度 、 重 度 灾害 发 生 频 
次 更 多 。 不 同时 期 干旱 灾害 重心 也 存在 变化 (图 
6) ,1958 一 1977 年 间 ,干旱 灾害 严重 区 分 布 在 西北 
部 的 张家口 市 .承德 市 及 北京 市 部 分 地 区 。1978 一 
1997 年 间 ,干旱 重心 主要 集中 在 南部 的 邯郸 市 、 邢 
台 市 .石家庄 市 .衡水 市 及 沧州 市 部 分 地 区 。 进 入 
21 世纪 以 后 干旱 重心 则 主要 在 中 部 的 保定 市 .衡水 
市 沧州 市 北京 市 唐山 市 .秦皇岛 市 及 天 津 市 等 部 
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图 4 京 津 避 地 区 SP16 变化 过 程 (1958 一 2017 年 ) 


Fig.4 Evolution of the SP/6 in the Beijing-Tianjin-Hebei region from 1958 to 2017 


图 5 


基于 Mann-Kendall 的 京 津 蔓 地 区 SPIO 突变 分 析 
Fig.5 Analysis of SPI6 abrupt change in the Beijing-Tianjin-Hebei region based on Mann-Kendall test 
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ABK, ENNA] BERETS h R 20 世纪 
90 年 代 ,干旱 重心 逐渐 向 东北 方向 跳跃 ,北部 沿海 
地 区 干旱 趋势 加 剧 。 这 与 本 文 结果 一 臻 ,干旱 重 心 呈 
现 出 由 北向 西南 ,在 向 中 东部 沿海 地 区 转移 的 趋势 。 
2.2 边缘 分 布 和 Copula 函数 选择 

干旱 历时 和 和 干旱 强度 的 边缘 分 布 拟 合 中 ,大 部 
分 站 点 最 优 分 布 为 指数 分 布 和 广义 极 值 分 布 ,少数 


(a) 1958 一 1977 年 


(b) 1978 一 1997 年 


km 中 0.06~0.15 


站 点 为 正 态 分 布 广义 帕 累 托 分 布 及 对 数 正 态 分 布 
等 。 经 验 分 布 和 理论 分 布 拟 合 效果 良好 (图 7) , 均 
能 通过 K-S 显著 性 水 平 检验 。 

基于 单 变量 边缘 分 布 , 构 建 Copula 联合 概率 分 
布 ,并 通过 AIC 检验 拟 合 优 度 ,确定 的 联合 分 布 函 数 
见 表 3( 部 分 站 点 )。 其 中 ,Frank Copula 拟 合 最 优 的 
站 点 31 个 , Gumbel Copula 拟 合 最 优 的 站 点 5 个 ， 


(c) 1998 一 2017 年 


图 例 
一 SPI6 等 值 线 
SPI6 


mm —0.25~-0.12 
mm —0.11~-0.06 
æm —0.05~-0.01 
æ 0~0.05 


国 0~0.06 
Lkm 0.07-0.17 


A6 RENAR 


ESEJE SP16 空间 分 布 


Fig.6 Spatial distribution of SP/6 in the Beijing-Tianjin-Hebei region with different periods 


10 - 
0.8 + 
, tee 
z = 
w 0.4 + K 
BK BR 
a — aRt 
K ---- 理论 分 布 
La 1 1 上 1 1 L 1 上 1 上 4 
0 123 465 © 7T 8 POT VB 0 2 4 6 8 10 12 14 16 18 20 22 
干旱 历时 /mon 干旱 强度 
到 7 干旱 历时 和 干旱 强度 累积 概率 分 布 ( 怀 来 站 ) 
Fig.7 Cumulative probability distribution of drought duration and drought severity ( Huailai station) 
表 3 Copula 函数 拟 合 优 度 检 验 
Tab.3 Goodness of fit test of Copula function 
ee y AIC AIC AIC 最 优 、 AIC AIC AIC 最 优 
F 站 名 (Clayton ) ( Frank ) ( Gumbel ) Copula F 站 名 ( Clayton) (Frank) (Gumbel) Copula 
1 张家口 一 198.81 一 279.33 -264.79 Frank 6 青龙 — 342.25 — 391.68 — 383.75 Frank 
2 HE -233.58 -271.20 -265.44 Frank 7 怀 来 — 258. 32 — 386. 26 -372.50 Frank 
3 承德 -281.92 -307.52 -280.31 Frank 8 定 州 -194.02 -264.60 -258.38 Frank 
4 霸 州 -329.19 -323.15 -298.54 Clayton 9 曹妃甸 -253.99 -301.38 -261.00 Frank 
5 盐 山 — 279.85 — 325.43 -303.09 Frank 10 深 州 — 293.56 — 363.81 — 323.65 Frank 
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Clayton Copula 拟 合 最 优 的 站 点 1 个。 结果 表明 , 京 
津 费 地 区 干旱 灾害 变量 的 联合 概率 拟 合 ,Frank Cop- 
ula 拟 合 效果 最 好 。 

2.3 ” 京 津 费 干 旱灾 害 单 指标 年 遇 型 分 布 

京 津 蔓 地 区 不 同年 遇 型 的 干旱 历时 .干旱 强度 
( 表 4, 部 分 站 点 ) ,表现 为 随 着 年 遇 型 增加 ,干旱 历 
时 和 干旱 强度 增 大 , 且 年 遇 型 越 大 ,指标 相差 越 明 
显 。5 a 一 遇 时 数据 差异 较 小 ,其 中 干旱 历时 最 小 值 
和 最 大 值 分 别 为 5. 31.7. 50 ,干旱 强度 最 小 和 最 大 
值 分 别 为 4.78 .8.41。 而 年 遇 型 变 大 ,数据 差异 变 
大 ,50 a 一 遇 干 旱 历 时 最 小 和 最 大 值 分 别 为 10. 82、 
21.16,50 a 一 遇 和 干旱 强度 最 小 和 最 大 值 分 别 为 
14.76 .83.24 ,差异 明显 。 

干旱 历时 空间 分 布 ( 图 8a ~8d) ,5 a 一 遇 历 时 
较 长 的 地 区 ,主要 分 布 在 中 部 和 南部 部 分 地 区 以 及 
东部 沿海 地 区 ,10 a .30 a、50 a 一 遇 历 时 较 长 的 地 区 
分 布 大 体 一 致 , 主要 位 于 京 津 囊 的 中 部 。 其 中 ,5 a 
一 遇 干 旱 历 时 较 长 地 区 为 保定 市 中 部 .廊坊 市 天津 
市 东部 .唐山 市 ,秦皇岛 市 、 承 德 市 南部 及 南部 的 邢 
台 市 中 部 . 邯 埋 市 中 部 地 区 。10 a、30 a 一 遇 时 , 干 
旱 历 时 较 长 地 区 为 保定 市 、 雄 安 新 区 EBT TH. 天津 
市 沿海 地 区 .唐山 市 南部 ,秦皇岛 市 沿海 地 区 及 张 家 
口 市 南部 等 部 分 地 区 。50 a 一 遇 历 时 较 长 地 区 为 保 
定 市 . 雄 安 新 区 张家口 市 中 南部 、 廊 坊 市 天津 市 南 
部 、 及 唐山 市 和 秦皇岛 市 部 分 地 区 等 。 

京 津 嘉 地 区 干旱 强度 分 布 (图 8e ~ 8h) 与 干旱 
历时 存在 差异 。 其 中 ,5 a 一 遇 水 平 下 ,承德 市 .唐山 
市 .秦皇岛 市 及 廊坊 市 存在 高 危险 区 ,10 a 一 遇 水 平 

下 ,张家口 市 西南 部 、 保 定 市 .廊坊 市 及 沧州 市 西北 
部 地 区 存在 高 危险 区 。30 a、50 a 一 遇 和 干旱 强度 较 


大 地 区 大 体 一 致 ,主要 位 于 中 西部 地 区 ,其 中 高 危险 
区 为 保定 市 西北 部 和 沧州 市 西北 部 地 区 。 
2.4 ” 京 津 费 干 旱灾 害 指标 联 合 分 布 

相同 年 遇 型 (5 a、10 a、30 a、50 a) 下 干旱 历时 
和 干旱 强度 的 联合 重 现 期 , 均 小 于 单 指标 ( 表 5) , 且 
在 50 a 一 遇 下 重 现 期 差异 较 大 ,空间 分 布 也 存在 差 
异 (图 9)。 对 于 联合 重 现 期 ,5 a、10 a 一 遇 水 平 下 ， 
邯郸 市 中 南部 .邢台 市 中 部 石家庄 市 东部 及 保定 市 
南部 .唐山 市 西部 .张家口 市 北部 地 区 ,联合 重 现 期 
较 小 ,为 高 危险 区 。 而 在 30 a、50 a 一 遇 水 平 下 H 
郸 市 .邢台 市 石家庄 市 .唐山 市 ,及 秦皇岛 市 . 张 家 
口 市 .保定 市 部 分 地 区 , 重 现 期 较 小 ,为 高 危险 区 ; 除 
了 南部 地 区 外 ,东部 沿海 及 张家口 市 西北 部 地 区 也 
易 发 生 严重 干旱 。 这 与 其 他 研究 结果 基本 一 致 , 即 
中 南部 为 干旱 多 发 区 '” ,呈现 南 多 北 少 , 蔓 南平 原 
最 多 '” | 。 同 现 重 现 期 表现 为 大 于 单 指标 重 现 期 , 且 
不 同 站 点 同 现 重 现 期 相差 较 大 ,高 危险 区 分 布 基 本 
一 致 。5 a、10 a 一 遇 水 平 下 ,张家口 市 西部 ,保定 市 
北部 廊坊 市 中 部 及 沿海 地 区 秦皇岛 市 和 沧州 市 , 同 
现 重 现 期 较 小 ,危险 性 高 。30 a、50 a 一 遇 水 平 下 ， 
表现 为 中 部 的 北京 市 沧州 市 .衡水 市 .张家口 市 西 
南部 及 承德 市 北部 地 区 , 同 现 重 现 期 较 小 ,为 高 危险 
区 。 

以 怀 来 站 联合 概率 等 值 线 为 例 (图 10a) ,等 值 
线 向 右 侧 弯曲 ,表明 干旱 历时 对 干旱 灾害 影响 较 大 。 
干旱 强度 为 0 ~5 时 ,等 值 线 分 布 密集 ,并 且 干 旱 历 
时 跨度 较 大 ,同时 干旱 历时 为 0 ~5 时 ,等 值 线 分 布 
也 较 密 集 ,主要 发 生 短 历 时 一 高 强度 .长 历时 一 低 强 
度 以 及 二 者 同步 的 干旱 灾害 。 此 外 ,将 不 同 干旱 灾 
害 类 型 的 实际 重 现 期 与 理论 重 现 期 进行 对 比 ,检验 


表 4 不 同 重 现 期 下 干旱 历时 和 干旱 强度 


Tab.4 Drought duration and drought severity of different return periods 


干旱 历时 干旱 强度 
序号 站 名 
T=5a T=10a T=30 a T=50 a T=5a T=10a T=30 a T=50 a 
1 张家口 6. 12 10. 08 16.36 19.27 6. 28 11.63 27.55 40. 36 
2 HE 6.07 10. 36 17.16 20. 32 7.64 13.04 21.59 25.57 
3 承德 6.04 9.67 15.42 18.09 8.41 11.91 15. 84 17.16 
4 TN 6.85 11.06 17.74 20. 85 7.45 12.03 19.29 22.67 
5 盐 山 5.39 8.13 14.02 17. 68 6.48 9.87 16.22 19.66 
6 青龙 7.50 9.11 11.06 11. 83 8.02 10.72 13.74 14.76 
7 怀 来 6.34 9.61 14.79 17.20 6.27 11.52 28.55 43.15 
8 EM 6. 84 11.15 17.98 21.16 7.37 10.05 14.30 16.37 
9 曹妃甸 6.53 10. 65 17.18 20.22 8.29 12.07 16.43 17.94 
10 PRI 5.63 7.99 12.51 15.04 7.15 10.18 14.64 16.61 
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(a) 5 a 一 遇 干 旱 历 时 (b) 10 a 一 遇 干 旱 历时 (c) 30 a 一 遇 干 旱 历 时 
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保定 ” 雄 安 新 区 
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(d) 50 a 一 遇 干 旱 历 时 (e) 5 a 一 遇 干 旱 强度 ®© 10 a 一 遇 干 旱 强度 
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沧州 


m 21.01~21.16 m 13.01~13.02 

= 18.01~21.00 = 12.01~13.00 

= 15.01~18.00 .01~8. = 11.01~12.00 

100 mm 12.01~15.00 0 100 = 6.01~7.00 0 Œ 10.01~11.00 
Lkm m10.84~12.00 L—km © 4.79~6.00 km © 9.32~10.00 


(g) 30 a 一 遇 干 旱 强 度 (h) 50 a 一 遇 干 旱 强度 
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一 等 值 线 一 等 值 线 
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图 8 不 同 重 现 期 下 干旱 历时 和 王 旱 强度 分 布 


Fig.8 Spatial distribution of drought duration and drought severity at different return periods 
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RS 干旱 特征 指标 联合 重 现 期 和 同 现 重 现 期 


Tab.5 Joint return periods and co-occurrence return periods of the drought characteristic index 


联合 重 现 期 同 现 重 现 期 
序号 站 名 

T=5a T=10a T=30 a T=50 a T=5a T=10a T=30 a T=50 a 
1 张家口 4.41 7.99 19.64 30.23 5.76 13.36 63. 45 144. 60 
2 HE 4.59 8. 50 21.08 32.10 5.49 12.15 52.03 113.05 
3 承德 4.37 7. 86 19.34 29. 84 5.85 13.74 66.85 154.01 
4 ai 4.37 7.60 18.48 28.73 5.85 14.61 79.73 192.39 
5 盐 山 4.33 7.77 19.12 29.58 5.92 14.03 69.55 161.46 
6 青龙 4.27 7.62 18.81 29.19 6.03 14.53 74.10 174.04 
7 怀 来 4.24 7.55 18.66 29.02 6.09 14.79 76.44 180.49 
8 定 州 4.17 7.40 18.36 28. 66 6.24 15.40 81.90 195.58 
9 曹 恕 馈 4.38 7.89 19.40 29. 93 5.83 13.65 66.09 151.89 
10 深 州 4.20 7.46 18.48 28.80 6.18 15.16 79.71 189.54 


联合 重 现 期 的 准确 性 。 以 河北 省 怀 来 站 为 例 ( 图 
10b) ,1958 一 2017 年 该 站 共 发 生 干 旱 事件 46 次 ,对 
这 些 干旱 事件 进行 分 类 吐 ] , 重 现 期 在 2 a 以 下 的 干 
早 事件 为 轻 度 干 旱 , 重 现 期 大 于 8 a 为 重度 干旱 ， 
2 ~8 a 则 为 中 度 干 旱 。 将 不 同 干旱 类 型 的 理论 重 现 
期 与 实际 重 现 期 进行 对 比 ( 表 6) ,结果 表明 绝 大 部 
分 在 灾害 事件 的 理论 重 现 期 与 实际 情况 符合 度 较 
高 ,具有 较 高 的 可 信和 度 。 实 际 重 现 期 与 理论 重 现 期 
的 Spearman 相关 系数 为 0.58( 通 过 95% 的 信和 度 检 
验 ) ,具有 良好 的 相关 性 ,说 明理 论 重 现 期 能 反映 京 
津 乙 地 区 实际 干旱 发 生 情 况 , 具 有 实际 意义 。 
结合 京 津 费 地 区 2015 年 粮食 单产 空间 分 布 
(图 11 ) 。 粮 食 单产 较 高 地 区 主要 分 布 在 京 津 费 的 
东南 部 ,包括 邯 韩 市 东部 .邢台 市 中 部 ` 石 家 庄市 东 
部 衡水 市 西部 .保定 市 东南 部 .北京 市 .唐山 市 及 承 
德 市 中 部 地 区 。 这 些 地 区 与 联合 重 现 期 高 危险 区 分 
布 基本 一 致 。 因 而 ,应 该 注意 干旱 灾害 预防 ,尤其 是 
春季 到 秋季 的 干旱 较 严 重 ,应 加 强 和 改进 河北 西南 
部 山区 农田 的 灌溉 系统 '”| ,整体 提高 干旱 防 控 水 


3 结论 与 讨论 


(1) 近 60 a 间 (1958 一 2017 年 ) 京 津 费 地 区 
SP16 变化 趋势 表现 为 ,20 世纪 50 年 代 末 上 升 趋势 ， 
20 世纪 60 年 代 初 下 降 趋势 ,1964 一 1966 年 又 转 为 
上 升 ,1966 年 至 今 整体 下 降 , 其 中 1979—1980 年 短 
暂 上 升 ;在 20 世纪 80 年 代 、1992 一 1994 年 间 、20 tH 
纪 90 年 代 末 ~21 世纪 前 10 a 下 降 趋势 最 显著 ; 存 
在 两 个 突变 时 期 ,分 别 出 现 在 20 世纪 50 年 代 末 ~ 


60 年 代 初 和 2015 年 附近 。 

(2) 干旱 历时 和 王 旱 强度 的 最 优 边缘 分 布 分 别 
为 指数 分 布 和 广义 极 值 分 布 ,其 中 Frank Copula 在 
计算 干旱 历时 和 干旱 强度 联合 分 布 时 拟 合 最 优 。 

(3) 京 津 束 地 区 干旱 灾害 历时 表现 为 中 西部 危 
险 性 较 高 ,主要 分 布 在 保定 市 廊坊 市 、 张 家 口 市 及 
北京 市 ,天津 市 .唐山 市 和 秦皇岛 市 也 有 部 分 地 区 危 
险 性 较 高 ;干旱 强度 表现 为 中 西部 地 区 危险 性 较 高 ， 
主要 分 布 在 保定 市 .张家口 市 沧州 市 及 北京 市 ,其 
中 雄 安 新 区 危险 性 也 较 高 。 

(4) 干旱 历时 和 王 旱 强度 的 联合 重 现 期 ,与 单 
旧 标 空间 分 布 存 在 差异 ;其 中 ,5 a、10 a 一 遇 下 , 唐 
山 市 石家庄 市、 保定 市 及 邢台 市 存在 高 危险 中 心 ; 
30 a\50 a 一 遇 下 ,北京 市 .沧州 市 .衡水 市 及 承德 市 
北部 等 地 区 的 重 现 期 长 ,危险 性 低 , 其 余地 区 较 短 ， 
危险 性 高 ; 同 现 重 现 期 在 5 a、10 a 一 遇 下 ,北京 市 、 
张家口 市 ,沧州 市 秦皇岛 市 以 及 承德 市 等 地 区 危险 
性 较 高 ,30 a、50 a 一 遇 下 , 则 表现 为 北京 市 .承德 市 
北部 沧州 市 以 及 衡水 市 等 地 区 危险 性 高 ,其 余地 区 
危险 性 低 。 此 外 ,站 点 的 理论 联合 重 现 期 和 实际 联 
合 重 现 期 较 符合 ,可 信和 度 高 。 


干旱 灾害 对 京 津 囊 地 区 的 粮食 生产 影响 较 明 
显 , 同 时 也 会 引发 饮水 困难 等 问题 。 了 解 干旱 灾害 


的 高 危险 区 以 及 干旱 重心 的 转移 变化 ,有 针对 性 地 
制定 不 同 致 灾 强 度 下 的 应 急 预 案 ,着 重 针 对 高 危险 
区 开展 防范 ,可 以 有 效 地 减轻 损失 。 未 来 干旱 灾害 
评 佑 研究 ,也 应 从 多 因素 ,多 维度 的 展开 ,使 评估 更 
具 科 学 性 。 

此 外 , 京 津 副 各 地 区 都 存在 不 同 程度 的 干旱 灾 
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Fig.9 Spatial distribution of drought duration and drought severity return period at different return periods 
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Fig. 10 Joint return period of drought duration and 
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表 6 不 同 干旱 历时 和 强度 的 理论 重 现 期 与 实际 
重 现 期 ( 怀 来 站 ) 
Tab.6 Theoretical and real return periods of different 


drought duration and severity (Huailai station) 
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1 1976.01 13 16.38 重 9.92 11.79 


2 1986. 06 6 5.91 轻 3.25 4.00 
3 1988. 06 2 2.19 轻 1.33 1.81 
4 1989. 12 7 5.85 中 4.75 4.26 
5 1997. 04 1 1.13 轻 0.83 1.38 
6 1997. 12 7 7.21 中 3.92 4.80 
7 1999. 03 3 3.12 轻 1. 92 2.27 
8 2010.01 8 6.59 中 3773 4.85 
9 2013.05 1 1.35 轻 0.58 1.44 
10 2014. 08 2 1.74 轻 1.25 1.62 
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Fig.11 Grain yield per unit area in the 


Beijing-Tianjin-Hebei region 


害 危 险 区 ,在 京 津 费 协同 发 展 背景 下 ,提升 干旱 灾害 
防 灾 减 灾 能 力 刻不容缓 。2018 年 4 月 4 河北 雄 安 新 
区 规划 纲要 》 中 指出 “加 强生 态 环境 建设 ,坚持 生态 
优先 .绿色 发 展 等 ”, 雄 安 新 区 及 其 周边 地 区 也 存在 
干旱 灾害 高 危险 区 。 这 对 干旱 灾害 预防 提出 了 新 的 
考验 ,针对 高 危险 区 ,不 仅 要 加 强 干 旱灾 害 管理 ,更 
要 注重 干旱 灾害 的 防治 ,这 样 才能 进一步 提升 京 津 
副 地 区 应 对 干旱 灾害 的 综合 能 力 。 
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Abstract; Drought disaster is one of the major meteorological disasters in the Beijing-Tianjin-Hebei region , China 
with serious economic and food security impacts. The hazard of drought disasters is evaluated to identify the high 
hazard areas,and provide references for the disaster prevention and emergency control. Based on the precipitation 
data of 37 meteorological stations from 1958 to 2017 ,the standardized precipitation index at 6-month scale ( SP/6) 
is calculated. And the drought duration and the drought severity are defined using the run theory. We use Kolmogor- 
ov-Smirnov method and A/C method to confirm the most fitted marginal distributions and the Copula function. The 
return periods of drought duration and drought severity are calculated under different return periods (RP 5 a,10 a, 
30 a,50 a). The results show that the SP/6 is in a downward trend from the mid-1960s to current times. The expo- 
nential distribution and the generalized extreme value distribution are the optimal marginal distributions of drought 
duration and drought intensity ,respectively. Frank Copula has the best fitting in the calculation of joint probability 
distribution. As far as the joint return periods of drought duration and drought severity show Beijing, Cangzhou, 
Hengshui and north Chengde of Hebei Province have a long joint return period with low hazard and the rest of the 
region has a short joint return period with high hazard ( RP 30 a,50 a). The theoretical joint return period of the 
drought disaster is close to the actual situation , which indicates that the joint return period can better reflect the situ- 
ations of drought disasters. 


Key words: drought; standardized precipitation index; Copula; joint distribution; Beijing-Tianjin-Hebei region 


